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ABSTRACT 

We present a new calculation of the intergalactic 7-ray pair-production absorption 
coefficient as a function of both energy and redshift up to the redshift of 3C279, z = 0.54. 
In reexamining this problem, we make use of new observational data on the intergalactic 
infrared radiation field (IIRF), together with recent theoretical models of the galactic 
spectral energy distributions of the IIRF from stars and dust reradiation and estimates 
of the IIRF from galaxy counts and COBE results. We present our results for two fairly 
well defined IIRF spectral energy distributions, one of which is within la of our previous 
estimate of the IIRF at ~ 20 ^m. We then apply our results to the 7-ray spectrum of Mrk 
421, and obtain good agreement with the observational data, including the recent results 
of the HEGRA group. 

Subject headings: 7-rays:theory - infrared:general - quasars:general - quasars: individual 
(Markarian 421, 3C279) 
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1. Introduction 

The magnitude and shape of the intergalactic infrared radiation field (IIRF) are fun- 
damental to understanding the early evolution of galaxies. The IIRF is determined by star 
formation processes occurring very early in the evolution of galaxies - emission of radia- 
tion from the photospheres of cool stars, dust reradiation, and redshifting to even longer 
wavelengths. The IIRF is especially sensitive to rates of early star formation, since much 
of the energy in starbursts escapes in the infrared. The study of the IIRF is therefore of 
great cosmological significance. 

We have previously pointed out (Stecker, De Jager and Salamon 1992 (Paper I)) that 
very high energy 7-ray beams from blazars can be used to measure the intergalactic infrared 
radiation field, since pair-production interactions of 7-rays with intergalactic IR photons 
will attenuate the high-energy ends of blazar spectra. Determining the intergalactic IR 
field, in turn, allows us to model the evolution of the galaxies which produce it. As energy 
thresholds are lowered in both existing and planned ground-based air Cherenkov light 
detectors (Cresti 1996), cutoffs in the 7-ray spectra of more distant blazars are expected, 
owing to extinction by the IIRF. These can be used to explore the redshift dependence of 
the IIRF. Furthermore, by using blazars for a determination of attenuation as a function 
of redshift, combined with a direct observation of the IR background from the DIRBE 
detector on COBE, one can, in principle, measure of the Hubble constant H at truly 
cosmological distances (Salamon, Stecker and De Jager 1994). 

The potential importance of the photon-photon pair-production process in high energy 
astrophysics was first pointed out by Nikishov (1961). Unfortunately, his early paper 
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overestimated the energy density of the IIRF by about three orders of magnitude. However, 
with the discovery of the cosmic microwave background radiation, Jelley (1966) and Gould 
and Schreder (1967) were quick to point out the opacity of the universe to photons of 
energy greater than 100 TeV. Stecker (1969) and Fazio and Stecker (1970) generalized 
these calculations to high redshifts, showing that photons originating at a redshift z will 
be absorbed above an energy of ~ 100(1 + z)~ 2 TeV. 

An EGRET-strength "grazar" (observed 7-ray blazar) with a hard spectrum extend- 
ing to multi-TeV energies is potentially detectable with ground-based telescopes (see the 
review on the Atmospheric Cherenkov Technique by Weekes (1988)). Pair- production in- 
teractions with the cosmic microwave background radiation (CBR) will not cut off the 
7-ray spectrum of a low-redshift source below an energy of ~ 100 TeV (see above). How- 
ever, as pointed out in Paper I, even bright blazars at moderate redshifts (z > 0.1) will 
suffer absorption at TeV energies owing to interactions with the IIRF. 

Shortly after Paper I was published, the Whipple Observatory group reported the 
discovery of ~ TeV 7-rays coming from the blazar Mrk 421, a BL Lac object with a 
redshift of 0.031 (Punch et al. 1992). This source has a hard, roughly E~ 2 , 7-ray spectrum 
extending to the TeV energy range (Mohanty, et al. 1993; Lin, et al. 1994) where the pair- 
production process involving infrared photons becomes relevant. Following this discovery, 
we made use of the TeV observations of Mrk 421 to carry out our suggestion of using such 
sources to probe the IIRF (Stecker and De Jager 1993 (Paper II); De Jager, Stecker and 
Salamon 1994). This problem was also discussed by Dwek and Slavin (1994) and Biller, et 
al. (1995). 
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There are now over 50 grazars which have been detected by the EGRET team (Thomp- 
son, et al. 1995). These sources, optically violent variable quasars and BL Lac objects, 
have been detected out to a redshift greater that 2. 

Of all of the blazars detected by EGRET, only the low-redshift BL Lac, Mrk 421, 
has been seen by the Whipple telescope. The fact that the Whipple team did not detect 
the much brighter EGRET source, 3C279, at TeV energies (Vacanti, et al. 1990, Kerrick, 
et al. 1993) is consistent with the predictions of a cutoff for a source at its much higher 
redshift of 0.54 (see Paper I). So too is the recent observation of two other very close BL 
Lacs (z < 0.05), viz., Mrk 501 (Quinn, et al. 1996) and 1ES2344+514 (Schubnell 1996) 
which were too faint at GeV energies to be seen by EGRET. 

In this paper, we calculate the absoption coefficient of intergalactic space using recent 
theoretical models of the spectral energy distribution (SED) of intergalactic low energy 
photons in conjunction with recent determinations of the IR background obtained by 
analysis of COBE data and galaxy counts. After giving our results on the 7-ray optical 
depth as a function of energy and redshift out to a redshift of 0.54, we apply our calculations 
to the high-energy 7-ray spectrum of Mrk 421 and compare our results with the published 
spectral data obtained by the HEGRA group (Petry, et al. 1996) and the preliminary 
spectrum obtained by the Whipple group (Mohanty,e£ al. 1993). 
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2. Absorption of High Energy Gamma- Rays by Low Energy Photons in Inter- 
galactic Space 

The formulae relevant to calculations involving the pair-production process are given 
and discussed in Paper I. For 7-rays in the TeV energy range, the pair-production cross 
section is maximized when the soft photon energy is in the infrared range: 

A(£ 7 ) ~ A e -^ = 2AE 1}TeV fxm (1) 

where A e = h/ (m e c) is the Compton wavelength of the electron. For a 1 TeV 7-ray, this cor- 
responds to a soft photon having a wavelength near the K-band (2.2 / um). (Pair-production 
interactions actually take place with photons over a range of wavelengths around the op- 
timal value as determined by the energy dependence of the cross section.) If the emission 
spectrum of an extragalactic source extends beyond 1 TeV, then the extragalactic infrared 
field should cut off the observed spectrum between ~ 10 GeV and a few TeV, depending 
on the redshift of the source. 

The TeV spectrum of Mrk 421 obtained by the Whipple Observatory group is con- 
sistent with an E~ 2 power-law, showing no significant absorption out to an energy of at 
least 3 TeV (Lin et al. 1994; Sreekumar, et al. 1996; Mohanty et al. 1993). It is this lack 
of absorption below 3 TeV which we (Stecker and De Jager 1993) used to put upper limits 
on the extragalactic infrared energy density in the ~ 1 fxm to 8 (im wavelength range (see 
also Dwek and Slavin 1994). 

Our upper limits on the IIRF, obtained from the Whipple data on Mrk 421 below 
3 TeV, ruled out various exotic mechanisms for producing larger fluxes, such as decaying 
particles, exploding stars, massive object and black hole models (Carr 1988). They are 
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consistent with the extragalactic near infrared background originating from ordinary stellar 
processes in galaxies (Stecker, Puget and Fazio 1977; Franceschini et al. 1994; see also the 
papers included in Dwek 1996). 

As pointed out by Dwek and Slavin (1994), the SED of the IIRF should consist of 
both starlight and dust reradiation components with the later component peaking in the far 
infrared region. Thus, the power-law SEDs which we used previously are approximations 
which are valid over a limited wavelength range. In this paper we utilize the SED models of 
Franceschini et al. 1994, together with recent data, in order to make a better determination 
of the IR background in intergalactic space and to calculate the 7-ray absorption coefficient 
from pair-production off the IIRF. 

3. The Intergalactic Infrared Radiation Field 

A theoretical upper limit to the intergalactic infrared energy density whose ultimate 
origin is stellar nucleosynthesis was given by Stecker et al. (1977). Direct observational 
limits on the IIRF, determined by analysis of the COBE/DIRBE data with other compo- 
nents subtracted out, has been given by Hauser (1996). Both of these determinations are 
well below the limits obtained by Biller, et al. (1995). Biller, et al. (1995) argued that the 
most conservative estimate of an upper limit on the IIRF would be obtained by taking 
the flattest spectral index for Mrk 421 consistent with EGRET data at GeV energies and 
assuming that the spectrum of Mohanty et al. (1993) is the result of absorption with 
an optical depth r > 6 in the TeV energy range. The huge IIRF required to obtain this 
optical depth is inconsistent with the COBE/DIRBE data (Hauser 1996). Also, in order to 
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produce a power-law absorbed 7-ray spectrum from a different power-law source spectrum 
in the TeV energy range, they had to employ an unphysical SED which is inconsistent with 
that produced by theoretical models of stellar and dust emission from galaxies. 

Recently, Puget, et al. (1996) have claimed a tentative detection of the far-infrared 
background from an analysis of COBE data at wavelengths greater than 200 ^m. An 
estimate of the IIRF from counts of IRAS galaxies has been given by Gregorich, et al. 
(1995). We have given both estimates of and upper limits to the IIRF (Stecker and De Jager 
1993; De Jager, Stecker and Salamon 1994) based on analyses of the preliminary spectrum 
of Mrk 421 at TeV energies obtained by the Whipple group (Mohanty, et al. 1993). All 
of these direct and indirect observational results can be brought together into a coherent 
scenario where the IIRF is assumed to be generated by stellar processes in galaxies. The 
resulting IIRF from both starlight and dust reradiation should then have a two-peaked 
SED made up of the contributions from all of the galaxies in the universe. The energy 
flux of the IIRF is expected to be within a factor of two of the value of 10 nWm _2 sr _1 . 
In Figure 1, we have plotted the observational data together with two theoretical SEDs 
of the galaxy-generated IIRF. These SEDs, the thick solid lines 1 and 2, are based on the 
models of Franceschini, et al. (1994). We used the starlight component exactly as given by 
these authors; However, the more uncertain dust component was normalized to a slightly 
higher value by using the data of Puget, et al. (1996) in the far-infrared range for their two 
different assumptions about subtraction of HII emission. The estimate of the IIRF from 
IRAS galaxy counts by Grigorich, et al. (1995) at 60 fim is consistent with both models 
within its uncertainty. 
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We also show an estimate of the SED at intermediate wavelengths from dust tori 
around Seyfert galaxy nuclei by Granato, et al. (1996). This additional component could 
add significantly to the lower SED model (model 1) at wavelengths between 10 //m and 
100 /im. 

Our chosen SEDs are generally within a factor of 2 of each other over the wavelength 
range of significance for the absorption of TeV 7-rays. We therefore feel that they represent, 
within the narrowest range of error presently available, a reasonable construction of the 
IIRF. We have used these two SEDs to calculate the opacity of intergalactic space to high 
energy 7-rays and to apply these calculations to observations of Mrk 421 in particular. 

4. The Opacity of Intergalactic Space to the IIRF 

The technique for calculating the pair-production optical depth as a function of energy 
has been given in Papers I and II. We do not here take account of changes in the IIRF with 
redshift, and instead make the assumption that the background is basically in place at a 
redshifts < 0.54, having been produced primarily at higher redshifts (Cowie, et al. 1996). 
We therefore limited our calculations to z < 0.54, the redshift of the powerful grazar 3C279 
which we discussed in Paper I. 

Figure 2 shows the results of our calculations of the optical depth for various energies 
and redshifts up to 0.54, the redshift of the powerful EGRET grazar 3C279. We take 
a Hubble constant of H a = 70 km s _1 Mpc _1 , which is consistent with various recent 
observations as presented at the recent Space Telescope Symposium on the Extragalactic 
Distance Scale. We assume for the IIRF, the two model SEDs given by the solid lines 
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marked 1 (lower curve) and 2 (upper curve) in Fig. 1. The two model SEDs yield different 
values for the optical depth r(E) only for E > 1 TeV, since the SEDs only differ at 
wavelengths greater than 2 //m. We note that r is greater than 1 for 3C279 for energies 
above ~100 GeV. Designing air Cherenkov telescopes with lower threshold energies should 
enable these sources to be detectable. 

We have obtained parametric expressions for Ti(E) and t 2 (E) (corresponding to mod- 
els 1 and 2 respectively) for z < 0.3. The double-peaked form of the SED of the IIRF 
requires a 4th order polynomial to reproduce parametrically. It is of the form 

4 

log 10 [Tj(E TeV , z)} = ^2{a i:j + b i:j log 10 z)(log 10 E TeV y for O.K E TeV < 50 (2) 

i=0 

for models j = 1 and 2. Table 1 gives the best-fit coefficients and bij for each model. 
These parametric expressions can be used by observers to correct for intergalactic absorp- 
tion and derive the true spectrum of sources out to a redshift of 0.3. 

Figure 3 shows spectra from various instruments during the quiescent phase of Mrk 421 
compared with the spectra calculated including absorption predictions for the two SEDs 
which we used for the IIRF. The EGRET differential spectrum of Mrk 421 (Sreekumar, 
et al. 1996) was converted to an integral spectrum in the range 70 MeV to 4 GeV, after 
adding the integral flux above 4 GeV to the EGRET differential points, assuming a model 
spectrum given by (for j = 1 and 2) 

J(E 7 ) = 1.3 x 10' 11 E-^expi-TjiE-^ = 0.031)) cm-V^eV" 1 . (3) 

Absorption does not influence the GeV part of the spectrum. The integral flux above 0.1 
GeV implied by eq. (3) is within 1 a of the published flux for Mrk 421 (Sreekumar, et al. 
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1996). By forcing this E~ x integral spectrum through the EGRET data, we obtain a value 
for x 2 of 7.5 for 5 degrees of freedom, showing consistency with our assumed spectrum. 
Normalizing our spectrum to the EGRET data in this way, we find excellent agreement 
with the TeV data as shown in Figure 3, even though the TeV spectra are from different 
time periods. 

Above 1 TeV, Petry, et al. have found a steep integral spectral index of 2.8 ± 0.6 for 
Mrk 421. The preliminary Whipple spectrum for Mrk 421 (Mohanty, et al. 1993) appears 
to be consistent with this steep index as shown in Figure 3, even though this spectrum is 
highly uncertain above 5 TeV (R.C. Lamb, private communication; Biller, et al. 1995). The 
spectrum given by eq. (3) connects the EGRET doXa. at GeV energies to the TeV spectrum 
of the HEGRA collaboration (Petry, et al. 1996) and follows the spectral steepening in the 
data above 1 TeV. 

The agreement between observation and calculation here supports our assumption 
that electron-photon cascading (Protheroe and Stanev 1993; Aharonian, Coppi and Volk 
1993) is not expected to significantly influence the TeV spectrum of Mrk 421. This is 
because the source is observed with a point-spread function of HWHM of ~ 0.2° and it 
can be shown that even an intergalactic magnetic field as small as 10~ 20 G would deflect 
the cascade electrons away from the source direction, creating a 7-ray halo much larger 
than the point-spread function. Since the measured magnetic field strengths in galaxy 
clusters and superclusters are of the order of a juG (Kronberg 1994 and refs. therein), it is 
highly unlikely that the intergalactic magnetic field would be as weak as 10~ 20 G. 
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5. Conclusions 

We have calculated the absorption coefficient of intergalactic space from interactions 
with low energy photons of the IIRF both as a function of energy and redshift, using 
reasonable theoretical SEDs which are consistent with present observational data. We 
have further applied our calculations to the grazar Mrk 421, which is the only extragalactic 
source for which spectral data presently exist at TeV energies. 

Observations of the 7-ray spectrum of the Crab Nebula from the same groups which 
have observed Mrk 421 both show agreement that the integral spectral index of the Crab 
Nebula is ~1.7 at ~ 1 TeV (Weekes, et al. 1994; Petry, et al. 1996). Furthermore, Petry, 
et al. (1996) have shown that the spectrum of the Crab Nebula above 1 TeV is flatter than 
the corresponding spectrum of Mrk 421. The index of the integral spectrum of Mrk 421 of 
2.8±0.6 derived by Petry, et al. (1996) is significantly steeper than the value of ~ 1 needed 
to connect the GeV and TeV data (see Figure 3) . While preliminary, the data of Mohanty, 
et al. (1993) are also consistent with this conclusion. 

Thus it appears that there is substantial observational evidence for a real steeping in 
the spectrum of Mrk 421 at a several TeV (see Figure 3) which may be naturally interpreted 
as an absorption cutoff caused by pair production interactions with the IIRF (De Jager, 
Stecker and Salamon 1994). For the two SEDs which we have chosen for the IIRF, the 
lower model (model 1) would allow for some additional intrinsic absorption in the source. 
However the use of either model will account for the spectral data (within the errors) 
without postulating any significant intrinsic absorption or a cutoff in the emission spectrum 
of the source at that particular energy. Both model SEDs give a value of r > 1 at 7-ray 

12 



energies above 5 TeV (see Figure 2). Therefore, we feel that absorption in intergalactic 
space by the IIRF is the simplest, most natural explanation for the TeV spectrum of Mrk 
421. It is an explanation which is consistent with present observations of the IIRF and the 
predicted absorption effect obtained using reasonable theoretical models of the IIRF. This, 
in turn, gives us confidence in our predictions of the absorption coefficient for sources at 
other redshifts. 

Our calculations confirm our previous conclusions in Paper I that the high energy 
(TeV) spectra of sources at redshifts higher than 0.1 should suffer significant absorption. 
The recent detections of two additional sources at redshifts below 0.1, viz., Mrk 501 (Quinn, 
et al. 1996) and 1ES2344+514 (Schnubell 1996), coupled with the non-observations of all 
of the other known EGRET blazars at redshifts above 0.1, even those which are much 
brighter than Mrk 421 in the GeV range and have comparable spectral indeces, further 
supports our calculations and interpretation of the Mrk 421 spectrum. 
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Table 1 



Coefficients for the polynomial in Eqn. (2) 
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-0.19 


0.07 



17 



Figures 



Fig.l 

The low energy intergalactic photon spectrum from the microwave to the optical, in- 
cluding the cosmic microwave background (CMB) spectrum and dust and starlight spectra. 
The theoretical SEDs given by the thick solid lines labled 1 and 2 are based on models 
of Franceschini, et al. (1994), slightly renormalized to the results of Puget, et al. 1996 
(filled circles). The filled square is the result of Grigorich, et al. (1995). The upper-limit 
arrows are the COBE/DIRBE residuals given by Hauser (1996), with the arrow lengths 
representing the range of the residuals. The lower-limit arrow near 2 fxm is from Tyson 
(1990), based on galaxy counts. We also show our previous upper limit (Stecker and De 
Jager 1993 - SD93) and estimate (De Jager, et al. 1994 - DSS94) based on the Mrk 421 
data from Mohanty (1993). The upper scale shows the characteristic energy for 7-ray 
absorption from photons of wavelngth A as given by eq. (1). 

Fig.2 

Optical depth versus energy for 7-rays originating at various redshifts obtained using 
the SED models 1 and 2 as shown in Fig. 1. 

Fig.3 

An E~ 2 power-law Mkr 421 spectrum with absorption calculated using SED models 
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1 and 2 of Figure 1 and fitted to the observational data from EGRET as described in the 
text. The Whipple data are from Mohanty, et al. (1993). The EGRET data are from 
Sreekumar, et al. (1996) and Sreekumar, private communication. The HEGRA data are 
from Petry, et al. (1996). The AIROBICC upper limit is from (Karle, et al. 1995), labled 
"A". 
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